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Abstract 



In this paper we study the asymptotic behavior for a nonlocal heat equation in an inho- 
mogenous medium: 

p{x)u t = J * u — u in x (0, oo) , 

where p is a continous positive function, u is nonnegative and J is a probability measure 
having finite second-order momentum. Depending on integrability conditions on the initial 
data Mo and p, we prove various isothermalisation results, i.e. u(t) converges to a constant 
state in the whole space. 



1 Introduction 

The aim of this paper is to study the asymptotic behavior for a nonlocal heat equation in an 
inhomogenous medium: 



p{x)u t = J * u — u, (x, t) € R N x (0, oo), 
u(x, 0) = uq(x), x € 



(1) 



Here, no is a nonnegative continuous function in Wl N and * denotes the convolution with a 
kernel J : R N -> R, which is a radial, continuous probability density having finite second-order 
momentum: 

I J(s)ds = l, E(J) = / sJ(s)ds = 0, V(J) = / s 2 J{s) ds < +oo . 
Jr n Jr n Jr n 

Typical examples of kernels that we consider are the gaussian law, the exponential law or any 
compactly supported kernel. We also assume that p is a positive, continuous function in 1^, 
whether integrable or not. 
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The operator J *u — u can be interpreted as a non-local diffusion operator. Indeed, if u(x, t) 
represents the density of a single poulation and J(x — y) is the probability to jump from y to 
x then the term (J * u)(x) is the rate at which individuals arrive to x and —u(x) is the rate at 
which individuals leave from x, see for instance [7]. In the case of heat propagation, u stands 
for a temperature and p(x) represents the density of the medium. 

Problem (pQ) is called non-local because the diffusion at u(x,t) depends on all the values of 
u in the support of J and not only of the value of u(x, t), as it is the case for the local diffusion 
problem 

f p(x)u t = Au, (x, t) E M. N x (0, oo), 
\ u(x,0) = u (x), x £ R N . 

For this local problem it is well know that for dimension N = 1, 2 there exists a unique solution 
in the class of bounded solutions, see [8] and [5]. Moreover, if p G L 1 (M Ar ) and uq is bounded, 
then as t — > oo the solution converges on compact sets to Ep(tto), the mean of uq with respect 
to p : 

/ My)p{y) d y 

JR N r- TO 



Ep(n ) := 



P(y) dy 



This phenomenon is called isothermalization, since the heat distribution converges to a non- 
trivial isothermal state in all the space. However, for dimension N 3 uniqueness is lost in the 
class of bounded solutions and some solutions decrease to zero as t — > oo, so that isothermaliza- 
tion does not take place, see [9] and [6]. 

On the contrary here, we show that if p is integrable, given any bounded initial data there 
exists a unique classical solution of problem (JTJ) and that the isothermalization effect always 
occurs for all dimension N ^ 1 (see more precise statements below). 

The case when p is not integrable is also considered, which is more related to the study of 
the homogeneous case (p = 1), see [H] and pQ. For bounded solutions, the flux at infinity is so 
big that solutions go down to zero asymptotically while if the data is unbounded, the solution 
may go to infinity asymptotically as t — > oo. 



Organisation and main results 

We first prove in Section [2] a comparison result which gives uniqueness for problem ([I]) . In 
Section Owe study the existence in the class of bounded solutions, which requires an approxima- 
tion first by functions p n that do not degenerate at infinity. The main theorem is the following: 



Theorem 1.1 Let p > 0, continous, and u$ be a bounded nonnegative continuous function. 
Then there exits a unique classical solution of problem ([T]). 

Before studying the asymptotic behaviour, an important step in this direction consists first in 
proving that if p is integrable and uo bounded, the following conservation law holds: 

/ p(x)u(x,t) dx = / p{x)uq{x) dx . 
Jr n Jr n 
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This is a nontrivial result which is obtained by approximation with Neuman problems in bounded 
domains, see Section UJ Then, Section [5] is devoted to study the isothermalisation phenomenon 
for bounded initial data. Then the main theorem is as follows: 

Theorem 1.2 Let p > 0, continuous, integrable and uq be a bounded nonnegative continuous 
function. Then u(x,t) — > E p («o) as t — > co in Li^ oc (M. N ) for any 1 ^ p < oo; the convergence 
also holds in L 1 ( y o): 

lim / \u(x, t) — Ep(uo)| p{x) dx = . 

Notice that we have no restriction on the space dimension contrary to what happens for the 
"local" heat equation. In the case when p is not integrable but still uq £ L x (p) (that is, 
puQ € L 1 (]R Ar )) the flux at infinity forces the solution to go to zero: 

Theorem 1.3 Let p > and uq be a bounded nonnegative continuous function such that uq E 
L 1 ( / o). // p is not integrable in R , then u(x, t) — > as t — > oo in \J^ oc {W N ) for any 1 ^ p < oo. 

Finally, in Section [6] we investigate the case when the initial data is not bounded. We first 
prove an existence result in the class of unbounded solutions provided p does not degenerate too 
rapidly at infinity. More precisely, if 

pW^TTW 7 ^ 2 ' (2) 

then we have an existence result for quadratic initial data: 

Theorem 1.4 Letuo be a positive continuous function with at most quadratic growth at infinity. 
If the function p satisfies ([2]) then there exists a minimal solution of problem ([T|). 

More generally, if p is integrable, we prove similar isothermalisation results for the minimal 
solution: 

Theorem 1.5 We assume that p is a continuous positive integrable function in M. N and that 
uq is a continuous nonnegative function, possibly unbounded, such that there exists a solution u. 
Noting u the minimal solution, the following holds: 

i) If u o £ L x (p) the isothermalization takes place in L 1 (/?) ; 

lim / \u(x, t) — Ep(uo)| p{x) dx = . 

ii) If uq (ji L 1 (p), we have that for all 1 ^ p < oo, 

lim u(x,t) = co in Lf (R N ). 

In the case of nonintegrable p's with uq £ L 1 (p) the asymptotic behavior is more difficult to 
treat. For instance, if p = 1, the solutions 

u(x, t) = \x\ 2 + V( J)t, and u(x,t) = 1 

have different behavior. Thus, there is a balance between p, J, and the initial data uq which is 
not easy to handle and the question remains open. 
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2 Preliminaries 



Let us specify first what is the notion of solution that we use: 

Definition 2.1 Let uq G Lj L oc (IR Ar ). By a stong solution of (|TJ) we mean a function u £ 
C°([0, oo); L 1 1 oc (IR'' v )) such that u t ,J *u € Lj^IR^ x (0,oo)), the equation is satisfied in the 
\j\ oc -sense and such that u(x,0) = u$(x) almost everywhere in M. N . 

We shall consider also solutions with more regularity: 

Definition 2.2 A classical solution of ([T|) is a solution such that moreover u,ut, J*u £ C°(1R' /V x 
[0, oo)) and the equation holds in the classical sense everywhere in R-^ x [0,oo). 

A classical sub or supersolution is defined as usual with inequalities instead of equalities in the 
equation. 

Now let us state a simple regularity Lemma, which contains a technical trick that we shall 
use several times in the sequel: 

Lemma 2.1 Let u be a strong solution of ([1]). We assume moreover that uo is continuous in 
M. N and that the convolution term J * u is continous in M. N x [0, oo). Then u and ut are also 
continuous in H N x [0, oo) and u is a classical solution. 



Proof. We introduce the following transform: 

T p [u](x,t) := e t/p{x) u{x,t). (3) 
A straightforward calculus show that v = T p [u] satisfies 

e */p(z) 

v t = — -r^— (J *u)(x,t), 
p{x) 

which is a continuous function in W N x [0, oo). Integrating between and t we get: 
v(x,t) = / dtv (x, s) ds + v(x, 0) = / dtv (x, s) ds + u(x, 0) , 



hence v is continuous in M. N x [0, oo). This implies that u is also continuous in M. N x [0, oo), and 
the equation holds in the classical sense. □ 



Remark 2.1 It is well-known in the convolution theory that under one of the following assump- 
tions, the convolution term is continuous: 

(i) u is bounded (since J is integrable); 

(ii) J compactly supported and u locally integrable. 



The following lemma concerns the comparison of classical sub/supersolutions of the problem. 
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Lemma 2.2 Let u and u be continous functions in R N x [0, oo). We assume that u is a classical 
supersolution of JU) and that u is a bounded classical subsolution of ([I]) with u(x,0) ^ u(x,0). 
Thenu^ u in R N x R, . 



Proof. We consider the functions 



S ,2 



ws = u- u - 5 - ^yTjs M , and <t>8 = T p [ws] ■ 



Since 



J*|x| 2 - \ x \ 2 = [ J(y)\y\ 2 dy-2 [ J(y) < x,y > dy = V( J) 

JR N JR N 

the above functions are related by the inequality 

p{x){4> s )t= e^W((J*^)(x,t) + |). (4) 

The function ws is continuous in R N x [0, oo) and notice that ws ^ —5 < at time t = 0. There 
are two options: 

0') either ^ -5/2 in 1^ x (0, oo); 

(ii) or there exists a point (x,t) € R^ x (0, oo) such that wg(x,i) > —5/2. 

Let us assume that we are in this second case and let to be defined as follows: 

t := inf{t > : 3x G R N , w s (x, t) > -5/2} < oo . 
A first estimate shows that 

ws(x,t)>-5/2^\x\<R s := (ili^NU) 17 



2 

< OO , 



so that the only possible points where ws may reach a level above —5/2 are located inside the 
fixed ball Br s . Thus, up to a finite number of terms, any minimizing sequence remains inside 
the compact set Br s x [0, to + 1] • After extraction and using the continuity of ws, we get that the 
inf is attained and there is a point (xo,to) such that w;^(xo,to) = —5/2. Of course this implies 
that to cannot be zero, since it; ,5 ^ —5 at t = 0. 

So, to > and necessarily ws ^ —5/2 in R N x (0,to), which implies also J * ws ^ —5/2 for 
t € (0, to), everywhere in R N . Using (jU), we see that (4>s)t ^ in (0,to) so that 

V(x, t)eR N x (0, t ) , 4>s(x, t) ^ <f> s (x, 0) = w s (x, 0)^-5. 

But taking x = xq and letting t — > to, by continuity we arrive at ws{xs,to) ^ 4>s(xo, to) ^ — 5 
which is a contradiction. 



So, we end up with the first possibility (i): for any 5 > 0, ws ^ —5/2 in all R x (0, 00) and 
taking limit as 5 — > we obtain 

(u-u)(x,t) for all x G R N , t > , 

which ends the proof. □ 
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Remark 2.2 Observe that the same result holds if we only assume that the subsolution grows 
strictly less than \x\ 2 : indeed, if u(x,t) ^ C(l + |x| 2 ~ £ ) for some e G (0,2) then the point xq 
remains inside a fixed ball with radius 0(5~~ 1 ^ £ ) and the rest of the proof follows identically. See 
also [lj. 

An immediate consequence of the comparison property is the uniqueness of the bounded 
classical solution. 

Corollary 2.1 Given uq continous and bounded in M^, there exists at most one classical solu- 
tion u of ([T]). 

Proof. Since uq is bounded, then the constant ||ito||<x> is a bounded (classical) supersolution. 
Hence, if there exists a (classical) solution u, using the comparison result we see that u is 
necessarily bounded by ||uo||oo- Then if we have two solutions u\ and v,2, since both are bounded 
we can again use the comparison result in both ways to get that u ^ v and v ^ u so that the 
solution is unique. □ 



3 Existence and uniqueness of bounded solutions 

We shall assume first that p is bounded below by some positive constant, so that we can apply 
the usual fixed-point approach to prove the existence of the solution for integrable initial data. 
Then we derive a general result for possibly degenerating p's by approximation. 

3.1 Non-degenerate case. 

Lemma 3.1 Let p ^ po > 0, continuous and uq G L 1 (M). Then there exists a unique strong 
solution u of © in the class C°([0, co); L 1 (R Ar )) . 

Proof. Following [3], we consider to > to be fixed later on and the space Xt Q = C°([0, to] ; L 1 (M Ar )) 
equipped with the norm = max {||to(-,t)|| L i ,0 < t < t }. Then we define an operator 
T : X to — ► X to as follows: 



J- Wo \ 



, (w)(x, t) := wq(x) H j {(J * w)(x, s) — w(x, s)\ ds 

Pn Jo > 

A straightforward calculus, using that p ^ po shows that: 

lT W0 (w) - T zo (z)\\ < \\w - z \\ L i + — - z\\ . 

Po 

We deduce that for any fixed initial data no G L 1 (]R Ar ), if to is sufficiently small, T uo is a 
contraction, hence there exists a unique solution u starting with u(0) = Uo, defined up to t = to: 

u(x, t) = uq(x) H — / {( J * u)(x, s) — u(x, s)\ ds , 
P Jo J 
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which implies that dtu exists as a h\ QC function, and that ([T]) holds in (0, io) i n the strong sense 
(not classical though, a priori). We then iterate the argument to construct a solution for all 
time t > 0. □ 



3.2 Degenerate case 

We construct here a solution of ([1]) for continuous and bounded initial data, using the approxi- 
mation p n — > p and also an approximation of initial data from L 1 (M Ar ). We consider the following 
problem: 

( p n (x)(u n ) t = J * u n - u n , (x,t) e R N x (0, oo), . . 

\ u n (x,0) = u (x)xn{x), x € R N , 

where by Xn we denote the indicator function of B n (the ball of radius n centered at the origin) 
and p n {x) = max{p(x); a n } where (a n ) is a sequence of positive real numbers strictly decreasing 
to zero. With this choice, p n ^ a n , p n decreases to p, and p n = p in any compact set for n 
sufficiently big (since p is continuous and positive). 

Since p n ^ a n > and itoXn £ L 1 (M Ar ) n L 1 (p n ) we can apply the results of the previous 
Subsection to obtain that there exist a classical solution u n . 

Now, in order to obtain a solution of problem ([1]) we pass to the limit as n — > oo. 

Proposition 3.1 Let p > 0, continuous anduo € L°°(lR Ar ) nC°(M Ar ) be a nonnegative function. 
For any positive decreasing sequence a n \ 0, let u n be the solution of ([5]) associated to p n = 
max{/9, a n }. Then the sequence (u n ) converges in L^ oc to the unique classical solution u E 
C ^ x [0,oo)) of flU) with initial data u . 

Proof. Since p n ^ a n > does not degenerate, we may use Lemma [37X1 to get a solution u n S 
C^fO,^)^ 1 ^)) with p n . We first observe that || uo||oo is a bounded classical supersolution 
of Q with ^ UQ n ^ 1 1 uo 1 1 oo so we may apply Lemma with p = p n to get that 

^ U n (x,t) ^ IKHoo • 

This allows to pass to the limit in L°°-weak* along a subsequence (still denoted by u n ): there 
exists a bounded function u such that u n u. 

We then deduce that J *u n — > J * u pointwise (since J is integrable) , and J * u is continuous 
since u is bounded (see Remark 12. ip . Moreover, J * u n is uniformly bounded so that the 
convergence also holds in L 1 1 oc (]R Ar x [0, oo)). On the other hand, p n dtu n converges weakly in the 
sense of distributions to pdtu, which would be enough to obtain a weak (distributional) solution. 
But we have more. 

Using transform T Pn introduced in Lemma |2"7T] we see that v n := T Pn [u n ] satisfies the equa- 
tion: 

e t/p n (x) 

d t V n (x,t) = r^( J * Un){x,t) , 

Pn{X) 
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hence dtv n converges pointwise and in Lj 1 (R x [0, oo)) to some / which is continuous. Inte- 
grating in time we get: 

v n (x,t) = Uon(x) + / d t v n (x,s)ds -> v(x,t) = Uq(x) + / f(x,s)ds, 
Jo Jo 

which is continuous in H N x [0, oo). 

Hence u n = e _t ' pn ^v n also converges pointwise and in Lj 1 (R^ x [0, oo)), and the limit is 
necessarily u, which is a continuous function in R^ x [0, oo). Coming back to the equation for 
u n , we see that every term converges strongly, and in the limit we have: 

put = J * u — u , 

where all the terms of the equation are in fact in C ^ x [0, oo)). Hence u is a classical solution 
of (H|) and Corollary 12.11 implies its uniqueness. The last step of the existence proof consists in 
saying that since any possible limit is in fact unique, then all the sequence converges to u, and 
there is no need to extract. □ 



Proof of Theorem II. 1L It is just the combination of Corollary 12.11 and Proposition 13.11 □ 



3.3 Approximation by Neumann problems 

Another way to get a solution of ([1]) in M. N consists in solving first the equation in the balls 
B n and pass to the limit as n — > oo. Since we are interested in getting a conservation law, it is 
natural to consider here the following Neuman problem, where Xn denotes the indicator of B n : 

p(x)d t u n = (u n (y,t) - u n (x,t)) J(x - y)dy , x G B n , 

U n (x,0) = U Xn(x) , 



Lemma 3.2 Let p > be a continuous function in W N and uq be a nonnegative continuous 
function such that uq € L 1 (p) n L°°. Then for any n ^ 0, there exists a unique solution 
u n £ C (0, T; L 1 (i? n )) of and u n is continuous. Moreover, the maximum principle holds 
which implies ^ u n ^ ||«o||oo an d the conservation law holds: 

I u n (x, t)p(x) dx = / uq{x)p{x) dx . 

J B n J B n 



Proof. Existence and uniqueness is done in [3] using a semi-group approach, in the case p = 1. 
Even if p is not constant here, it is bounded from below in B n so that the modifications are similar 
to those of the non-degenerate case, section I3TT1 Moreover, using the same techniques as in [2], 
the equation (or inequation for sub/super solutions) still holds at the boundary, dB n x (0, oo). 
The comparison principle is obtained in the same way as in lemma I2T21 without the term |x| 2 . 
Therefore, using the constant solutions c = and c = ||«o||oo to compare with the solution u n , 
we get the desired estimate. Finally, the conservation law is simply obtained by integrating the 
equation over B n x [0, t\. □ 
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Proposition 3.2 Let p > be a continuous function inM. anduQ be a nonnegative continuous 
function such that uq € L 1 ( y o) nL°° and for any integer n > 0, let u n be the solution constructed 
in Lemma \3.Sh Then as n — > oo ; along a subsequence u n — > u in Lj 0C (M JV x [0,oo)) where u is 
the unique classical solution of ([T]). 

Proof. Since the sequence u n is bounded, it converges (along a subsequence still denoted u n ) 
in L°°-weak* to some nonnegative and bounded function it. Since J £ L 1 , this implies that 
J * u n — > J * u strongly and that dfU n — > dtu in the sense of distributions. We can then pass to 
the limit in the sense of distributions but here also we want a better convergence. 

Let us introduce a modified version of transform T as follows: 

v n (x,t) := e^ J * Xn ^P^w(x,t). 

Since ([6]) can be written as 

p(x)d t (u n )(x, t) = [J * (u n Xn)] (x, t) - (J * Xn){x)u n {x) , 

it follows immediately that v n satisfies the equation 

d t V n = p-r J * U n . 

p{x) 

This implies that dtv n converges strongly on compact sets of ~R. N x [0, oo), and so does 
v n {x,t) = uo(x)xn{x) + /o ( v n)t(x, s) ds. Then u n also converges strongly on compact sets of 
R N x [0, oo) to its limit u. 

Passing to the limit in the equation, we see that u is a strong solution of ([T]), which implies 
that it is the unique classical solution of this equation - this follows from Lemma 12.11 and 
Corollary 12. li Any other converging subsequence leads to the same solution u so that all the 
sequence u n converges to u in L 1 1 oc (R Ar x [0, oo)). □ 

Remark 3.1 By uniqueness of classical solutions, we have that the classical solution obtained 
by approximation of the function p{x) and the classical solution obtained as the limit of solutions 
of Neumann problems are identical. 



4 Conservation Law 

In this Section we investigate various conditions under which solutions of ([1]) satisfy a conser- 
vation law. A first lemma in this direction is as follows: 

Lemma 4.1 Let p be a nonnegative measurable function and u a strong solution of ([1]) such 
that u G C°([0, oo); L 1 (M Ar )) . If u (x) € L x (p) then for any t > 0, u(-,t) £ L 1 ( y o) and 

/ p(x)u(x, t) dx = / p{x)uq{x) dx. 
Js. N Jr n 
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Proof. First, we integrate equation ([!]) with respect to the time variable to obtain 



p(x)u(x, t) = p(x)uq{x) + I (J * u — u)(x,t) dt 



o 

We integrate over Br using Fubini's theorem: 

p(x)u(x, t) dx + I I u(x, t) dx dt = / p(x)uo(x) dx + / / J * u(x, t) dx dt . 
b r Jo Jb r Jb r Jo Jb r 

Using the fact that uq(x) G L 1 (p) and u G C°([0, oo); L 1 (R Ar )) , we can use monotone convergence 
for the right-hand side terms which implies that u(-,t) G L 1 (p). Moreover, since u(-,t) G L 1 , 



(J * u)(x,t) dx = / u(x,t)dx, for any t > 
Jrn 

p(x)u(x,t) dx = / p(x)uq(z 
Jr n 



so that finally we get 

(x) dx . 

□ 

In particular, this result is valid for the solutions constructed in Section 13. 11 but it is difficult 
in general to be sure that u(t) remains in L 1 . So we use the approximation by Neuman problems 
in order the keep the conservation law in the limit when p G L 1 . If p is not integrable, we have 
only a one-sided estimate but this will be sufficient to obtain the asymptotic behaviour. 

Proposition 4.1 Let p > 0, continuous and u G C ^) n L 00 ^) n L x (p). Then 
(i) for any t > 0, u(-,t) G L 1 ( y o) and 



p(x)u(x, t) dx ^ / p(x)uo(x)dx. 
Jrn 

(ii) if moreover we assume that p G L 1 (]R Ar ) then the conservation law holds: 

p(x)u(x,t) dx = / p(x)uo(x) dx . 



Proof. Since u G L°° n C°, we may use the approximating sequence {u n } of Section \3.3\ which 
satisfies that u n — > u at least almost everywhere, with ^ u n ^ \\uo\\ and 

u n (x,t)xn(x)p(x)dx = u (x)xn(x)p(x)dx . (7) 

Since no G L 1 (p), the dominated convergence theorem yields the convergence of the right-hand 
side integral as n — > oo. For the left hand side, using Fatou's Lemma we obtain that 

p(x)u(x, t) dx ^ liminf / p(x)u n (x,t)x n {x) dx = / p(x)uo(x) dx , 

rwo ° Jm Jr n 

which proves assumption (i) of the proposition. 

Finally, if we assume that p G L 1 ^) we can use also the dominated convergence theorem 
for the sequence u n XnP, which is bounded by ||ito||ooP £ L ■ We then pass to the limit at the 
left-hand side of and get (ii). □ 
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5 Asymptotic behaviour for bounded solutions 

We shall now derive our main results concerning the asymptotic behaviour for ([1]). We divide 
the proof in several steps. 

5.1 Weak limit 

This first step is easy, it only comes from the fact that the solution is globally bounded: 

Lemma 5.1 Let p > 0, continuous anduo G C (]R Ar )nL oo (R Ar ). Let u be the associated classical 
solution. Then for any s > there exists a subsequence tk — ► +oo such that the following limit 
exists in L°° -weak*: 

u 00 (x,s):= lim u{x, s + ifc). 

t k ^oo 

Proof. Since no is bounded, u is also bounded by comparison in the class of classical solutions, 
thus there exists a subsequence tk — > oo such that u(-, s+t^) converges in L°°-weak* to a function 

« 00 (-,s)eL 00 (l JV ). □ 

5.2 Lyapounov functional 

We now want a stronger result, so we use a Lyapounov functional: 

Lemma 5.2 Assume the hypotheses of Lemma \5.1\ and that uo G h 2 (p). Then there exists a 
constant C = C(uq,p) such that, 

poo p 

p(x)(ut) 2 (x, s) ds ^ C. 

Proof. First we prove that for the approximating problem ([5]) the following functional 
F[u n ](t) = / J(x - y)(u n (x,t) - u n (y,t)) 2 dxdy 



is a Lyapunov functional which is nonincreasing along the evolution orbits. Multiply equation 
([5]) by (u n )t and integrating in space we obtain that 

4-F[un](t) = -4 f Pn (x)((u n ) t ) 2 (x,t)dx, (8) 
dt J r n 

while multiplying the equation by u n gives 

F[u n ](t) = -2^- [ p n (x)u 2 n (x,t)dx. (9) 
at j RN 

Notice that derivation under the integral is possible at the level n since we know that u n (t) G 
C^O.oojL 1 ^)), which implies that also d t u n G C(0,oo;L 1 ). 
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Integrating ([9|) we have that for some C = C'(uq, p), 

F[u n ](s) ds = 2 / p n (x)ul(x)xn{x) dx - 2 / p n (x)u 2 n (x, t) dx < C' . 
Indeed, as uo € L 2 (p) we have by monotone convergence that 

PM )X^^I dx<oo. 



Hence, i h- > i ? [n n ](t) is a decreasing function which is in 1/(0, £) for all t > 0. It follows that 
F[n n ](t) must be bounded for all f > by some constant C = C(uo,p). 

Moreover, F[u n ](t) is positive so that, integrating ([8]), we get for any t > 0: 

p((u n ) t ) 2 (x,s)dxds < -.F[u n ](t) ^ C(uo,p). 

Using Fatou's Lemma and the fact that p n (u n )t converges strongly to put, we obtain the desired 
result. □ 

As and immediate consequence of this result we obtain: 

Lemma 5.3 Assume the hypotheses of Lemma \5.1\ and that uq £ L 2 (p). For all sequence 
tk — ► oo and s > 0, 

|| \fp{^) u(-,s + tk) - \/pD w(-,tfe)||L2( R jv) -> as n -> oo . 

Hence, the limit function u^^Xjs) does noi depend on the variable s > 0. 

Proof. Note that for all sequence tk — ► oo, we get 

||v / pO)^(-,s + *fe) - V / P( : ) n (-^fc)||L2(RiV) = / /" «t(a;,o-)dcr) dx 



^s/ / p(x)(ut) 2 (x, a) da dx — ► 0. 



□ 



5.3 The cu— limit set 

We define the limit set as follows 

lj(uo) = {uoo G C°(M. N ) : 3tj — > oo such that u(-,tj) — ► Uoo(-) in L°°-weak*} 

Lemma 5.4 Under the hypotheses of Lemma \5.1\ and that uq £ L 2 (/?), i/ie u -limit set is reduced 
to constants. 
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Proof. Since u(x, s + tk) converges weakly in L°°-weak*, then as tk — > oo, 

(J*n)(x,s + t fc ) = / J(x-y)u(y, s + t k ) dy pm -^ lsc (J *u 00 )(x,s) . 



Moreover, since u is bounded, the convergence of J *u is also strong in Lj . On the other hand, 
dtpu(s + tk) — > pd s u ao (s) in the sense of distributions. We then pass to the limit in the sense of 
distributions in the equation and get 

d 

p(x)—U 00 (x, S) = J * Uoo(x, s) - Uoo{x, S) . 
OS 

Using Lemma 15.31 we know that Mqo is independent of s so that bounded solution (in 

the sense of distributions) of 

J * u oo — u oo = in R . 



We deduce that u^ is continuous because the convolution term is continous, and so it is nec- 
essarily constant in all Mr. Indeed, the equation for u^ implies that it is harmonic (see for 
instance H] for a proof) and bounded harmonic functions are constant. □ 

5.4 Identification of the limit 

We are now ready to identify the w-limit set. 

Lemma 5.5 We assume the hypotheses of Lemma \5.1\ and that uq € L 1 (p). Then the following 
holds: 

(i) ifpeL 1 (R N ),Lo(u ) = {E p (u )}; 

(ii) tfp$L 1 (m N ),Lo(u ) = {0}. 

Proof. Notice first that since uq € L 1 (p) n L°°, then uq € L 2 (p), hence we may use Lemma 
15.41 In the integrable case, p £ L , we observe that as u is uniformly bounded, the dominated 
convergence Theorem gives 



p(x)u n (x, s + tj) dx — > / p{ 



x) dx . 



Therefore, by Proposition 14. II - (ii) we obtain u^, = Ep(tio), so that the w-limit set is reduced to 
{Ep(n )}. 

In the case p L 1 (]R Ar ), we take a compact set K such that 

p(x) dx > / p(x)uq(x) dx , 



K JB. N 

which is always possible since uq £ L x (p). Using Proposition 14. H - (z) . Lemma 15.41 and Fatou's 
Lemma, we obtain 

^oo \k p( x ) & x ^ liminf \ K p(x)u n (x, s + tj) dx 

n^oo 

^ liminf f wN p(x)u n (x, s + tj) dx 
< f RN p(x)u (x) dx . 
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This implies that necessarily Uoo = 0, hence the w-limit set is reduced to {O}. □ 
5.5 Proofs of Theorems D and Q 

As consequence of the fact that the u— limit is given by only one function we can pass to the 
limit in the time variable without extracting any subsequence. Then it only remains to check 
that the convergence is better, which is done by using transform T p . 

Proof of Theorem 11.31 Under the hypotheses of the Theorem, we have that uq is continuous, 
bounded, and p not integrable, but nevertheless uq € L x (p) nL 2 (p). Thus Lemmas 15.11 and 15.51 
imply that for any s > 0, at least along a subsequence t n — > oo we have u(x,s + t n ) — > in 
L°°-weak*. But the same arguments are valid for any other subsequence such that u(x, s + t' n ) 
converges weakly. Since the limit is always zero, we deduce that for any s > 0, 

L°°-wcak* 

u(x, s + t) ► , 

t— »oo 

which implies that (J * u(s + i)) converges strongly in as t —* oo. Then, 

L°°-weak* 

p(x)d s u(x, s + t) = (J * u(s + t)) (x) — u(x, s + t) > . 

Even more, from lemma l5\2l we obtain that for any compact set K, 

/o(a;)[itt|(a;,s)da;^ ds < (J p(x)\u t \ 2 (x, s) dx^j (J p{x) dx^j ds < C(K, u Q , p). 

Then, at least for some sequence i& — > +oo, we have p(x)d s u(x, s + tk) — > in Lj 1 ^. Summing 
up, we obtain that 

lim u(x, s + tk) = in Lj 1 ^ . 

t fc ->+oo 

Of course, if t t— ► u(x,s + t) were to converge in L loc along another subsequence — > oo, the 
limit would necessarily be zero, so that finally u(-, t) — > in L 1 1 oc (M Ar ) as t — > oo. Moreover, since 
t i— ► u(-,t) remains bounded in L°°(]R Ar ), we deduce that the convergence holds in L^ )C (lR Ar ) for 
any 1 ^ p < oo. □ 

Proof of Theorem 11.21 The first part is done exactly as in the proof of Theorem 11.31 except 
that p is integrable here so that the limit is not zero, but M p (uo). To end the proof in this case, 
it only remains to prove the L 1 (/o) convergence. We fix e > and choose R > big enough so 
that (remember that p is integrable): 



J\x\>R 



(x) dx ^ e . 

<\x\>R 



Then 



\u(x, t) — u 00 \p{x) dx ^ 2e||iio||oo + / \u(x,t) — u oc \p(x)dx , 
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and using the L* convergence we get: 

limsup / \u(x, t) — u 00 \p(x) dx ^ 2e||uo||oo • 

t^oo JM.N 

Since e is arbitrary, we get that the limit is zero, which ends the proof. □ 

6 Unbounded solutions 

In this section we derive some results for unbounded initial data and solutions. Let us mention 
that in the case p = 1, further results are to be found in [I]. But here we still face the problem 
of the space inhomogenity implied by p. 

6.1 An existence result for quadratic growth 



We assume that there exists some rj > and ^ 7 ^ 2 such that for any x 6 



p(x) ^ — -\— . (10) 

This hypothesis means that p does not degenarate too rapidly at infinity. We then begin by 
constructing a supersolution of the problem (recall that V( J) is the second-order momentum of 
J, which is a finite constant): 

Lemma 6.1 // (fTQ|) holds, then for any A ^ Y(J)/r/ and A > the following function is a 
(classical) supersolution of ([T]): 

U(x,t) := Ae xt (l + \x\ 2 ). (11) 

Proof. It is a simple calculation: dtU = \Ae xt {\ + |x 2 |), 

J*U-U= Ae xt (J* \x\ 2 - \x\ 2 ) 

= Ae xt (j RN J(y) |y| 2 dy-2 j RN J(y) <x,y> dy) 
= Ae xt V(J). 

Hence we have a supersolution provided Xp(x)(l + |x| 7 ) ^ V(J). Since by assumption we have 
p{x){l + \x\ 1 ) ^ rj, it is enough to impose Xij ^ V(J), hence the result. □ 

Proposition 6.1 Let us assume that ()10p holds. Then for any nonnegative uq € C°(lR Ar ) ; 
satisfying uq(x) ^ C(l + \x\ 2 ) for some C > 0, there exist a strong solution u of JT]) with 
u(x, 0) = uq(x). 



Proof. Let us first consider an approximation uq u = uq ■ Xn where Xn is smooth, nonnegative, 
compactly supported and Xn /* 1- Let u n be the unique solution of ([T]) with initial data uq u 
given by Proposition 13. 1^ then by applying the comparison result for bounded solutions, the 
sequence u n is nondecr easing. 
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On the other hand, if A is big enough, and for instance A = V(J)/ry, we may use the 
supersolution U defined in (fTTj) to compare with u n . Notice that U is not bounded, but this is 
allowed in Lemma 12.21 which gives: 



(x,t) ^ Ae xt/v {l + \x\ 2 ). 



Hence the sequence u n converges to some u and we are able to pass to the limit in J * u n by 
dominated convergence, using that U is integrable with respect to translations of J. 

Using now Lemma [2TT1 we deduce that u is a classical solution of (pQ) and the initial data of 



Remark 6.1 This construction does in fact give a minimal solution: if u\ is any other solution, 
then it can be used as a supersolution for any u n and passing to the limit shows that u ^ u±. 
One can think that if we restrict the initial data to grow at most like |x| 2 ~ e , then uniqueness 
holds because the comparison argument is valid in this class (see Remark I2.2() . However, it is not 
clear whether the constructed solution enters this class unless we know more about uo, see [T]. 

6.2 Asymptotic behaviour for unbounded solutions when p is integrable 

We prove now that if u$ is integrable with respect to p, the isothermalisation phenomenon occurs 
(whether infinite or not). Notice that we gave sufficient conditions for existence of a minimal 
solution in the previous section. The first result is the following: 

Proposition 6.2 Let u$ E C°(M Ar ) (~lL 1 (p), p € L 1 (R Ar ) and assume there exists a solution u 
such that u(x,0) = uq(x). Then, if u denotes the minimal solution, we have 



Proof. Let u be the minimal solution and let us use the same monotone approximations that 
were used in Proposition 16.11 Since u n is bounded, Lemma 12.21 implies that we have a bound 
from above: 



Moreover, the convergence of u n to u is monotone, so we can pass to the limit in the above 
equation to obtain 



u is Uq. 



□ 




But since u n (x,0) £ 



u(x, t) ^ u n (x, t) in R x [0, oo) . 
L°°(R Ar ), Proposition O implies 





Using the above three equations we get that 
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where, 

\u(x,t) — u n (x,t)\p(x) dx = I {u(x, 0) — u n (x, 0))p(x) dx . 



\u n (x,t) - E p (uoXn)\p(x) dx, 

h = I \^ P (uoXn) -E p (u )\p(x)dx. 
Jm. n 

Observe that I\ and I3 are independents of t and tend to zero as n — > 00. Moreover, u n satisfies 
the hypothesis of Theorem 11.21 so that I2 tends to zero as t — ► 00. Therefore, we first have 



limsup / \u(x,t) — E p (uo)|p(x) dx ^ I\ + I3 
Jrn 



so that taking the limit as n — > 00 yields the desired result. □ 
In the case when uq £ L 1 (p), then infinite isothermalisation occurs: 



Proposition 6.3 Let p G L 1 ^) and u G C ^) such that u g L x (p). Then for any 
solution u with initial data uq and any 1 ^ p < 00, the following asymptotic behaviour holds: 

lim u(x,t) = +00 in L?(R N ) . 

t— >+oo 



Proof. As before, if there exists a solution, then we can approximate the minimal solution u 
by the family u n used in Proposition 16. 1[ Since this approximation is monotone, 

u(x, t) ^ u n (x, t) in R N x [0, 00) . 

But u n (x,0) satisfies the hypotheses of Theorem 11.21 so that 

liminf u(x, t) ^ lim u n (x,t)=c n , 

where c n = K p (u n (x,0)), the limit holding in all L^ )c (M Ar ). Hence passing to the limit as 
n — > +00, we obtain the result for u since c n — > ¥, p (uo) = +00, thus the same holds for any 
other solution. □ 

Theorems 11.41 and 11.51 follow from the conjunction of Propositions and 16. 1\ 16.21 and 16.31 
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